ITEP-TH-40/01
O
O
(N
<
oo
in
d
RELATIVISTIC TODA CHAIN AT ROOT OF UNITY III.
RELATIVISTIC TODA CHAIN HIERARCHY
S. PAKULIAK AND S. SERGEEV
Abstract. The hierarchy of the classical nonlinear integrable equations associated with
^ ' relativistic Toda chain model is considered. It is formulated for the iV-powcrs of the quantum
operators of the corresponding quantum integrable models. Following the ideas of the
paper pi[ it is shown how one can obtain such a system from 2D Toda lattice system.
The reduction procedure is described explicitly. The soliton solutions for the relativistic
Toda chain are constructed using results of |I2] in terms of the rational tau- functions. The
vanishing properties of these tau-functions are investigated.
^ . Introduction
en
O ' This paper is a supplement to the papers pl, 0] . It contains proofs of several statements
O ■ used in these papers and related to the classical relativistic Toda chain (RTCh) model. The
O \ papers [jl], 0] deal with the quantum integrable model based on the Weyl algebra at the roots
of unity
a'
X.
H , Namely, the simplest possible model which is the RTCh model at the root of unity was
considered there. The parameter lo being different from unity describes the deviation of the
quantum model from the corresponding classical counterpart.
In the quantum integrable models of this type, in contrast to the models where the quan-
tization parameter is in general position, there exists a non-standard classical limit. In usual,
in order to get a classical model from the quantum one, the quantization parameter is sent
Key words and phrases, relativistic Toda chain, classical integrable models, bilinear formalism.
This work was supported in part by the grant INTAS OPEN 00-00055. S.P.'s work was supported by
the grant RFBR 01-01-00539 and grant for support of scientific schools RFBR 00-15-96557, S.S's work was
supported by the grant RFBR 01-01-00201.
2 S. PAKULIAK AND S. SERGEEV
to unity and the Poisson bracket appears in this classical limit as normalized commutator.
In this limit the quantum observables becomes classical ones, namely, C-valued functions
depending on the evolution parameters (times or coordinates). In the quantum integrable
models based on the Weyl algebra at the root of unity there is another possibility. One may
consider the A^-th powers of the quantum dynamical variables
(2)
N N
Writ
ing the deformation parameter in the form
(3)
^ ^ ^2ni/N_^e/N2
and considering the limit e — > |^ one may obtain from the commutation relations (|I|) the
Poisson bracket
(4) {Um, Wn} = SnmUmWn ■
It is clear that logarithms of the dynamical variables form the canonical symplectic form
(5) n = y^ d \n Um /\ d \n Wm ■
The important feature of this unusual classical limit is that the quantum variables do not
disappear and their A^-th powers being the parameters of the unitary representations of the
Weyl algebra form the classical integrable system. The quantum analysis of such a combined
systems becomes more transparent. In particular, in the paper |^ a modified Q-operator
was constructed for the quantum RTCh model at the roots of unity using several result from
the corresponding classical integrable system.
The aim of this paper is to prove some of these results from the point of view of the
standard theory of integrable nonlinear differential equations. In this approach the integrable
equation appears usually as a representative of some family of equations sometimes called
hierarchy of integrable equations. In some cases, these hierarchies are certain reductions of
a most general system of equations. The subject of the present paper is the RTCh equation
and the associated hierarchy. It is known that this model can be obtained as a reduction of
the 2D Toda lattice (2DTL) system |jTT], |T2[ and its integrable discrete analogs
Relativistic generalizations of the classical integrable equations attracted much attention
after seminal work by S. Ruijsenaars [in] and have been investigated in several papers (see, for
3
example, 0, |^, |^ ) • The bilinear formalism for one of most simple example of these relativistic
generalizations, namely, for RTCh equation was presented in [^]. In the latter paper the
M-soliton solution to this equation were constructed in terms of Casorati determinants.
We show that these solutions coincide identically with the standard soliton tau-functions
constructed in KP or 2DTL theories [^, |16|, |T^ providing that parameters of solitons lies
on the rational curve, which is related to Baxter's curve appeared in the quantum integrable
models based on the Weyl algebra at the root of unity.
The paper is organized as follows. In Section |] we consider RTCh model in finite volume
(with lattice size equal to M) and using L-operator approach construct RTCh hierarchy in
the framework of slightly modified AKNS construction. In Section S we obtain the RTCh
hierarchy as certain reduction from 2DTL theory following ideas of the paper |TT[ . In last
Section the soliton solutions to the first equation of RTCh hierarchy are considered and iden-
tified with those used in papers ||1|, 0. The Appendix contains a short review of Casoratian
technique, which helped authors of the paper |]12l to find soliton solutions of RTCh equation.
1. RTCh hierarchy in t-matrix formalism
In order to establish relation with the papers [H, 0] we construct in this Section the
RTCh hierarchy using language of L-operators and r-matrix formalism. We will show that
RTCh hierarchy can be obtained as slightly modified AKNS system. In particular, we will
explain the meaning of the constant shift in the corresponding trigonometric classical r-
matrix mentioned in the paper [||.
The L-operator for RTCh in the form of 2 x 2 matrix
(6) Lm
y/UmWm \ ZWm
appeared first in |p. In (|^) z is a spectral parameter, k G C is the complex parameter and Um,
Wmi m E Tj satisfy the local Poisson bracket (^). This L-operatorQ is associated with lattice
site m. In this section we will consider the lattice of the finite volume, so m = 1, . . . , M.
^We used here a gauge-equivalent dependence of the L-operator ^ on the spectral parameter in contrast
to [|l], H and also changed its normalization.
4 S. PAKULIAK AND S. SERGEEV
It is easy to show p| that the Poisson bracket (H) can be rewritten in r-matrix form
(7)
{Lm{z) f Ln{y)} = Snm.[r{z, y) , Ln{z) ® L„Xy)]
where r{z, y) = r{z, y) + s and
r{z,y)
f a{z,y) \
(3{z,y)
(3{z,y)
V aiz,y) )
1 z^
y^
P{z.y)
zy
is the standard trigonometric classical r-matrix for £[2 and
/ \
-1/2
(9)
1/2
/
is a constant shift. One may show that after quantization this constant shift is related
to twisting of XXZ quantum i?-matrix in the commutation relation of the corresponding
quantum L-operators (see @, 0). In the quantum case the twisting of the quantum R-
matrix is necessary in order to use minimal finite-dimensional representations of the Weyl
algebra at the root of unity. On the classical level such a shift of the classical r-matrix
will simplify the corresponding classical integrable hierarchy by removing higher non-linear
terms (see discussion around (p^).
In order to introduce the evolution of L-operator (^) with respect to continuous times one
have to consider the monodromy matrix
(10)
Li L2
^M
A{z^) zB{z^)
zC{z^) D{z^)
where v4(2^), B{z'^), Clz"^) and D{z'^) are polynomials with respect to z'^ of the orders M,
M— 1, M—1 and M—2 respectively. The commutation relation (|^ shows that the coefficients
M
k=0
Tfc, k = 0, . . . , M of the trace of the inonodromy matrix ( |I0| )
M
(11) T{z') = A{z') + D{z') = J2Tk
are in involution with respect to this Poisson bracket
{Tfc,T„} =
and can be taken as Hamiltonians which correspond to evolution of the L-operator (^ with
respect to times t^
(12) ^ = {L^,T,}.
Our goal is to rewrite the equation ([12|) in the zero curvature form. To do this we will
use a method described in the book ||^ . We have
(13) {Lra{z),T{y)} = Vm-iiz,y)Lmiz) - L^{z)Vm{z,y) ,
where
Vm{z,y) = tT2[{id® Li{y))-- -{id® Lm{y))f{z,y) x
(14) X (id ® Lrr,+i{y)) ■ ■ ■ (id ® LM{y))]
and tr2{A ^ B) = tr(i?) ■ A. Moreover, since the singular part of the r- matrix (^ is propor-
tional to the permutation operator the following equality holds
lim (2:^ - y^)[Vm-i{z,y) Lm{z) - Lm{z)Vm{z,y)] =
which shows that r.h.s. of (|13|) is polynomial with respect to z and y
M
(15) VUz,y) = Y.Qli^\z)y'\
k=0
where
k
(16) Q^r^\z) = Y,<^,,{T)Q^^\z),
p=0
S. PAKULIAK AND S. SERGEEV
, , / ^2" /2
\ -z2"/2
(17) V
^ / /l2n/2
\ -/l2„/2 ^
and the entries of the matrix polynomials Qm (z) h2k, ^2k-i and f2k-i are differential poly-
nomials of the functions u^ and Wm+i only. The coefficients a^piT) in decomposition (jT^)
are polynomial functions of the integrals of motions T^.
We can see that the structure of 2 x 2 matrix polynomials ([17D looks very similar to those
in the AKNS construction. Let us introduce new time variables xj- such that evolution of
L-operator (|) is given by the zero curvature equation
(18) ^^ + Q^ZUz) LUz) - LUz) Q^:i\z) = .
In (^) we can choose integer n to be any non-negative integer and call the obtained
hierarchy of equations as RTCh hierarchy. Let us write explicitly several first equations of
this hierarchy (see the Section y, formula ( piP , where linear combination of the systems ([T
and (|20|) is identified with RTCh equation):
- dlnum _ dlnwm _ ^
dxQ dxo
(20)
OX I
(21)
dX2
d\\iw.
m
dxo
dxi
= UmWmWm+l{KUm + Mm-l) +
+ Um^iWm-lWm{K.Um-l + Um~2) + wl^^KUm + Mm-l)^ ,
The compatibility condition between the systems (^) and ( pTj ) falls into quite simple form:
(
dUm d'^Urr, „ OUr
'"m „ ^ ""m
f22)
0X2 OXl OX I
0X2 oxf Ox^
7
As we already mentioned, the structure of the matrices Qm (0) shows that the whole
RTCh hierarchy can be obtained as a version of AKNS hierarchy |^, |Tl| with small modi-
fication related to the shift matrix s ^. In (0) we decompose the matrices Qm into two
pieces. The first one correspond to the r-matrix (^ and to standard AKNS construction,
while the second one to the constant shift matrix s (|^). If the parameter e in (^) is free then
the first equation in the corresponding AKNS construction, written as the zero- curvature
condition for matrices (p!7|),
has the form (we renamed Um = e and Wm+i = / for simplicity)
(24)
* - S%4„/|i + 2(2.+ l) e^|^ + (.+ l)e3f
dx2 dxf dx^ \ dx^
Sf - «V+4«/|^ + 2(2.H-l)f/'|l-(. + l)eV'
dx2 dx\ dx^ \ dx^
It is seen from (0) that these nonlinear equations are simplified essentially at the choice
of the parameter e = — 1/2 by removing the higher nonlinear terms. The same phenomena
occurs with all higher hierarchy equations. We would like to repeat once more that such a
simplification is a classical analog of the possibility to use the minimal representation of the
Weyl algebra in the corresponding quantum model.
Due to locality property {Qm depends only on Um and Wm+i) the AKNS construction
do not describe discrete time evolution (plsl) . As well as in the usual AKNS hierarchy the
dependence on the discrete time or the discrete evolution (|l^) can be reconstructed after
using the Backlund-Schlezinger transformation [l^. The L-operator (|]) plays the role of
matrix which generates this transformation. It is possible to develop this AKNS construction
in full details using the methods described in the book 0, but instead we will describe in
the next section the RTCh hierarchy as a reduction of 2DTL.
2. RTCh as reduction of 2DTL
In this section we recall the basic facts on 2DTL model following the fundamental paper
T^. Then we consider a certain reduction of this system which leads to the RTCh hierarchy.
S. PAKULIAK AND S. SERGEEV
2.1. 2D Toda lattice. In |1J] 2DTL was defined in terms of Z x Z matrices L and L of
the following form
(25) L = ^diag[6j(m)]A^' , C = J^ diag[cj(m)]A^' ,
where bi{m) = 1, C-i{m) ^ for any m G Z and Am,k = Sm+i,k is a shift matrix. In
components, these infinite matrices can be written as follows
oo
^m,k = Sm+l,k + / ^b^j{m)Sm-j,k ,
(26) _ ^=° oo
Cm,k = C-i{m)Sm-l,k + ^Cj{m)Sm+j,k ■
For arbitrary infinite matrix A = ^.ggdiag[aj(m)]A'' we define its projections on the
upper-diagonal (including diagonal) A^ and the strictly lower- diagonal (excluding diagonal)
A_ matrices by the formulas
(27) A+ = Y^ diag[aj(m)]A-'' , A_ =Y^ diag[aj(m)]A^' .
i>o i<o
2DTL is a system of differential equations for the functions b-j{m), Cj{m), j > and
c_i(m) with respect to two copies of the continuous parameters (times) x = {xi,X2, . . . }
and y = {yi, 2/2, • • • } of the form
f28)
where Bg = {'C'^)_^_ and C^ = (£ )_.
The system ( P^ can be obtained as the compatibility conditions of the spectral problems
(29) CW{x,y) = W{x,y)A, ZW{x,y) =W{x,y)A-^
and the linear problems
dy^Wix,y) = C,W(x,y), dy^W{x,y) = CsW{x,y) ,
9
where W(a;, y) and W{x, y) are Baker- Akhiezer functions. Writing them in the form
^{x,y) = J^diag[i/'j(m;a;,2/)]A ^ ■exp(^(a;,A))
i=o
>V(x,y) = ^diag[^j(m;x,?/)]A^ ■exp(^(y,A ^))
j=0
W{x,y)-' = exp(-e(x,A))-^A--'diag[^;(m + l;x,y)]
3=0
oo
yV{x,y)-^ = exp(-e(i/,A ^))-^A^diag[^*(m + l;a;,i/)],
i=o
with iljQ{m;x,y) = ipQ{m;x,y) = 1 and ipQ{m;x,y) ^ for m G Z and using a bihnear
formahsm one can prove |114[ that the scalar generating series
-^{m;x,y;X) = A'"exp^(x, A) ^ ^'^(m; x,?/)A ^ ,
j=0
oo
^*(m;x,?/;A) = A~™exp^(-x, A) ^ V'*(m; a;, y)A^
(31) _ L='_
■^{m;x,y;X) = X"" exp^{y, \~^)^'iljj{m;x,y)X^ .
j=0
oo
-^*{m;x,y;X) = X ™exp^(-y, A ^) ^ V^*(m; x,y)A^
i=o
can be expressed in terms of family of tau- function r„(x, y), n E 1^
<i/{m;x,y;X) = ^ ,
^*im;x,y;X) = I^A^lI^^ ^
(32) >^'''^\^^^
Tt7^ \^ 'T^m+i[x,y-e[X))
^{m-x,y-X) = ^ ,
"^ {m;x,y;X) = ^ .
rm{x,y)
In ( |3lD the function ^(x, A) is
oo
e(x,A) = 5^x,A^
s=l
and in ( ^2]) the notation £:(A) means the set {A, A^/2, A^/3, . . . }.
The formulas ( |5^ and the spectral problems ( POD allow to express all the entries of the
matrices C and C as certain combinations of tau-function Tn{x,y). In particular, the first
10 S. PAKULIAK AND S. SERGEEV
coefficient c_i(n) in the decomposition of the matrix L has the form
(33)
c_i(m;a;,2/)
2.2. RTCh hierarchy as reduction from 2DTL. Consider the following substitution for
the entries of the matrices L and L (j > 0) W% '-
(34)
b^j{m)
Cj[7n)
c_ilmj
[K,
3 3
j=l j=0
+ u.
m—j—l) 1
"^n^'^+^n^'n^
m+iK'^'^m+j ' '^m+j+l) '
i=0
i=0
""m—1
Ur.
In ( |3^ ) we introduced the functions Um{x,y) and Wm{x,y) which depend on the continuous
times Xs and ys and k G C is a parameter, the same as in the previous Section.
Comparing the last formula in (^^ with ( p3D we can choose the following parametrization
of the function Um{x, y) in terms of tau-function
'Tmi.X.y)
(35)
Umix,y)
r.
m+l
'x,y)
We consider the similar ansatz for the variable Wm{x,y)
(36)
Wm{x,y)
e.
m+l
[x,y)
Om{x,y)
The formulas (|35|), (|3^) and (|3^) allows to get simple expressions relations for the matrices
C and £ in terms of RTCh dynamical variables Um and Wm- We have
Proposition 1.
(37)
C
e
r
w
r
l + A-
1
-U{A-k)
1-K-^A
where U , W, T and are diagonal matrices
U-' (1+A"M W
T
6 '
-l^ TJ/-1
T
e '
(38) f/^i = diagf-u
.±11
W = diag[u;,
±11
T
±1 ^ diag[r;;^i] , Q^' = diag[C]
and t/ie rational functions
1 — X
in formulas (W^ is always understood as the series
^ oo
— — = Vx^'
1- X ^-^
j=0
11
Corollary 1. The matrices C and C which correspond to RTCh reduction satisfy the relation
(39)
c-c = c-c
-K
Proof. Substituting formulas ( ^4|) into (126|) we obtain
l-A-'C = (I + k)
(40)
'CA + K= = (1 + k)
where we introduced the diagonal matrices
e
r_ 1 + A-i T^
'r K- A r '
(41)
r± = diag[r„±i] , 6 = diag[6'„] , 6*^
Tm+lTm-l + I'^'^n
(1 + K)em
Excluding from formulas (^0[) the diagonal matrix B we can find first the relation (|39|) and
using this relation in (|40|) we prove the statement of the Proposition. D
At the first sight the formulas (pTj) does not depend just on variables m„ and Wn, since
they includes explicitly the diagonal matrices T and 0. But this is not true since, we have
obvious formulas
(42)
e
A^'
T
r
H
e
A^'
T
r
H
diag
diag
n
.4 =
n
i=l
'^m+i^m+i
'^m—i'^m—i
j>0
j<0
Using formulas ( pTf ) we can obtain explicit formulas for the elements of the matrices Bs
and Cs which describe the evolutions of the matrices C and C by the formulas (pH]). This
will give the evolution equations of the functions Um{x,y) and Wm{x,y) with respect to
continuous times Xg and ys- Form of these equations is given by the following
Proposition 2. The substitution (^[j describes the reduction of 2DTL to the couple of the
RTCh hierarchies with respect to the sets of the times Xg and yg- These hierarchies are given
by the formulas
dlnum
(43)
dxg
d In Wr,
dxg
E IfiJ.. '"+'-E
3r
Jk>0 i=0
Jk>0 i=0
r=l
m — I
+E>
r=l
12 S. PAKULIAK AND S. SERGEEV
and
r_^s9lnuj
s-1
dys
(44)
where
(45)
and
dhiw
\s m
dys
i^.>o j=o \ r:
jfc>o j=0
Uj(m) = (-)-'if;^ (k(1 - (5j,o)um-i+i + ^im-i)
(46)
Wj(m) = (-)%„^ {k{1 - (5j-,o)w-i^.+i + w-i^.) .
The meaning of the Proposition ^ can be understood as follows. Take the equations from
the system (|28|) , which describe the evolution of the operator C and C with respect to the
time Xs- This set contains infinite number of equations, which describe the evolution of the
entries b_j {m) and Cj {m) with respect to this time. The statement of the Proposition |^ means
that the whole set of these equations reduces after substitution ( p^D just to pair of equations
given by the formulas (|^). The same is true for the evolutions with respect to times ys- One
may verify that first two pair of equations from ( |13| ) which correspond to s = 1,2 coincide
precisely with equations (|20| ) and (^) obtained from L-operator approach. A general proof
that the hierarchy (|T8|) obtained in L-operator approach coincide with hierarchy (^) for any
s should be a repetition of the analogous proof for the equivalence of such a presentations
in case of standard AKNS hierarchy. We do not go into these details here.
Explicit formulas (|43|) and (^ij) yield the following
Corollary 2. The systems (^) and j^^) transforms to each other after the simultaneous
transformations
(47) Xs < — > -{-K^ys , Un < — > w~^ .
The first equation in j^^ ) goes into the second equation of (^^ while the second of 1 ^^ ) to
first of g^.
13
The transformation (^) allows to forget about dependence of our dynamical variables on
times ys and consider only the set of equations (^31).
This symmetry of the reduced 2DTL is related to certain Backlund transformations of
2DTL investigated in [|T7| and used in the paper []T^ in order to obtain the bilinear formula-
tion of RTCh equation. Ended, formulas (^0|) contains 'new' tau-function dm given by (l4ll)
as a certain rational combination of r^ and Om- It is actually unclear from our presenta-
tion whether the functions dm and 6m are also the tau-f unctions. The last equation of the
substitution (0) and (|33D allows to identify only r^ with some tau-function of 2DTL.
Nevertheless, this assertion is true and one can use the following argument to see that the
function 9m is also a tau-function. It is clear that the symmetry transformation p7|) moves
the substitution fBll) into another one:
3 J
c_i(m) = ,
Wm-l
where now the ratio of functions Wm should be identified with a ratio of tau-functions ac-
cording to the formulas (^). The ansatz ( |5BD allows to identify the function 6m with some
another tau-function of 2DTL. This tau-function should be different since otherwise the
substitutions (^) and (^Sj) becomes trivial (the product UmWm = 1 for any m G Z in this
case).
On the other hand if one uses the substitution (|48|) in order to get expressions for the
matrices C and C similar to (^OD a new rational combination of tau-functions
^ ' "^ " (1 + ^)rm
(k)
will appear. Let us denote all these tau-functions as a set Tm
We can extend the value of the index k to be any integer number by repeating the transfor-
mations given by the formulas ( Pj) and ppj). It was shown in [jT^ that the whole set of these
tau-functions satisfy a reduced discrete 2DTL (the transformations (ETl) and (|i9|) describe
14 S. PAKULIAK AND S. SERGEEV
the evolution with respect to discrete parameter k) and all of them are some tau-functions
of original 2DTL hierarchy.
(k)
We will see in the next Section that in the soliton sectors the tau-functions Tm differs
by the proportional change of the soliton amplitudes (see (p3|)). In terms of the set Tm the
transformation (147|) can be written in the form
(50) r!^)^r!^-^)^ keZ
3. Soliton solutions of the RTCh equation
m
In this Section we will give the solitonic solution of the RTCh equation obtained in [12
terms of rational tau-functions. We will identify them with tau-functions used in the papers
|jl|, 0] as well as with soliton tau-functions appeared in KP |15, 16| and in 2DTL [14] theories.
In order to do this we have to transform the RTCh equation (120) into bilinear form analogous
to one used in the [jT2[. The authors of this paper used so called Casoratian technique to
find a solitonic solutions to the RTCh equation in form of the Casorati determinants. For
the readers convenience we will describe part of this technique relevant for our paper in the
Appendix.
We start with rewriting the equations (|20|) into bilinear form. Actually, in order to use
results on the reduction described in the previous section, we have to consider certain linear
combinations of equations from the RTCh hierarchy. Instead of (pOD we consider a linear
combination of (0) and (pOD
■9 In Mm , \ f , ,\
— = Wm{K,Um + Um-l) - [H + I) ,
(51) .1
7{ = Um{HWm + Wm+lj - [K + I) .
OX I
such that the system (pT|) has a constant {um = w^ = const) solution.
Excluding from (^) Wm and replacing Um = e'^"^''^^^^^^^'^^ we obtain (g^ = dqm/dxi
.52^ " = (Vr + ^ + g.:„_i)(v/r + ^ + gJ _ {VI + k + q'J{Vl + K + q'^^^)
The latter is the relativistic Toda chain equation. In the non-relativistic limit k — > cxd we
obtain from (^) the standard Toda chain equation
„" plm-l—gm. „qm — qm + l
Hm ^ c
15
Proposition 3. |]T2[ The substitution
'^m\-^lj
transform the equations ( |31| ^ into bilinear form
Dl{Tm O Trn) = 2(1 + k) (6',„+i^m_i - T^) ,
(54) -Dl(rm O Tm-l) = (1 + K){TmTm-l - OmOm-l) ,
(i + KjUm.t'm ^ Tm+lTm—l + ^''"m •
In ( |S^) the Hirota operator D{f o g) = f'g — fg' have been used.
We can also find a bihnear relation between tau-functions r^ and 6*^
(55) -01(6*^-1 O Tm) = (1 + K){0mrm-1 - ^Im-irm)
which is a consequence of the RTCh equation and Backlund transformation r^ -^ 9^ dis-
cussed at the end of the previous section. In the Appendix we will show how this bilinear
equation can be solved using Casoratian technique.
Let A* and A are two independent complex parameters which belong to the rational curve
(56) K = A*A-A*-A.
In the paper [|l| we used the following uniformization of the curve ( ^61)
A^ = e"^'^ ( a/cos2 (/) + k + cos
(57) ^
A^ = e*"^ f ycos^0 + K + cos 4,
which is convenient for restricting the model to the chain of the finite size with periodic
boundary conditions. For real k and the parameters A* and A are complex conjugated.
Define functions
(i-A:r(i-A^-i)V^-^p"i
^ ' (i-a^)'"(i-a;~i)Vi-^^)"i
f (1 + /t)W2
16 S. PAKULIAK AND S. SERGEEV
which depend on discrete time m, first continuous time Xi and some additional integer
parameter k. Define also the M-soliton tau-function by the determinant^
^!n^ (xi, {0fc,/fc}|fli) = det ip^^^^_{{xi;(j)jjj)
\ / «j=i
One of the statements of the paper [|T2[ is the following
(59)
Proposition 4. |T2| The tau-functions which solve the equations ( [5^j are given by the fol-
lowing specialization of the function^^)
,M
(60)
Tm [Xi; {(t)k, fk}\k=l) = Tm [Xi] {(pk, fk}\tLij ,
Om [Xi; {(j)k, fk}\tLl) = Tm^^ [Xi; {(l)k, fk}\kLl) ,
Om(xi;{^kJk}\k=i) = ^mMxi; {0fc,/fc}|Jlj .
The tau-functions given by the Proposition || coincide with those used in the papers |jl|, Q
up to certain normalization factor. Let us change this normalization in order to fix precise
relation to formulas of these papers and also to demonstrate the connection of these tau-
functions with standard solitonic tau-functions known from the 2DTL theory [|14|. We obtain
(A* = A* and A, = A<^J
(61)
where
(62)
and
(63)
-A,r-
r^' = Il(^: - a;) n Ti^v ^" ^^^-^' ^^L^^
2<J
TL^^ = ^(^*-^*)"'det
i<j
A* M
i,j=^
X^^\x,;<j>J) = fe^^*^-^^^^^
A*
m+l
w)r , <</>)
a; (1 - A
A
1-a;)
It is clear that factor before Tm in formula (^) cannot spoil the form of the equations (54)
so can be dropped. Now the tau-functions Tm = 'Tm and 6^ = Tm coincide identically
with those used in papers [||, If. The only difference is the shift in amplitude (|63|), but this
is because we used Um = Tm/Tm+i instead of Um = Tm-i/Tm in |l], H-
^The superscript '^^^ counts the number of Backlund transformations (p5|), not the number of sohtons as
in the papers (ij ||.
17
The formula ( |62D may be rewritten in the form analogous to the standard presentation of
solitonic tau-functions appeared in KP ||T3| or 2DTL theories [|T1| with all higher continuous
times but Xi are freezed:
M
A * A
(64) T^;^^ = det
6ij - X^^ {xi ; (f)j ,fj)Yl ^hk
M. ^'-^^
*J = 1
where
(65) s,,u = s{ct>A'), K^,^0 = |^ ^!_^f
In order to come from formula ( ]62| ) to (^4|) one have to use explicitly the inverse Vandermonde
matrix and this is a reason of getting product of factors Sj^k in soliton amplitudes. Actually,
the determinant presentation (p^ ) is equivalent to the recurrent relation which relates the
tau-functions with different number of solitons:
M-l
fc=i
where the phase shift is given by cross-ratio
.... , _ (a:-a-)(a.-a,)
^ ' " (A*-A,)(A,-A*) •
Specializing the M-soliton tau-function Tm {xi, {(j)k, fk}\kLi) to the zero value of contin-
uous time Xi = and setting the amplitudes /^ = 1 we obtain that it becomes proportional
to
(67) T^^^O;{^kA}\ti) - det\i-y-' {A:^^--' - A,^^--%^^ ■
This proves the following
Proposition 5. The rational solitonic tau function Tm (0; {0^, 1}|^{^]^) vanishes for values
of the discrete time m = —M + 1, —M + 2, ... ,0 irrespective the values of the spectral
parameters (pk-
This property of the rational tau-function was crucial in [Q in order to construct modified
Q-operators.
18 s. pakuliak and s. sergeev
4. Discussion
So far the discrete variable m in the previous sections runs over all integer numbers. It
is clear how to apply the results of the previous section to finite volume system or how to
satisfy the periodic boundary condition with respect to discrete time m
(68) Um[Xn) = Um+MyXji) , Wm[Xn) = Wm+M[Xn) ■
One should specialize the spectral variable which uniformize the curve (|56| ) to the finite
set
(69) '^" = 1^^ fc = l,...,M-l.
The case k = or k = M are excluded from (p9| ) because at this value of the parameter 0,
A = A^ and the solitons disappear. It is clear that in this case there is only finite number
(M including a zero soliton case) rational tau-functions.
Acknowledgment
The authors are grateful to G. von Gehlen, S. Kharchev, E. Sklyanin and F. Smirnov for
useful discussions and comments.
S.P. would also like to thank Max-Planck Instiut fiir Mathematik (Bonn) for support and
hospitality and S.S. would thank the hospitality of MPIM during his short visit to Bonn
supported by Heisenberg-Landau program.
Appendix: Casoratian technique
The goal of this Appendix is to demonstrate Casoratian technique |T^ and prove the
bilinear relation (^) with tau-function given by the Casorati determinant (|59|). Let us fix
values of the discrete parameters m and u. The functions ( |55D satisfy the following dispersion
relations
r% (m+i/)
(70) dx, ""^ "'^ ^"'+''+' '
V i + « [Vm+k - '^m+k ) - ^m+k+l ' ^ ^ /, .
im+v)
19
Let (^^'' means the column
(71)
' (Pm^''\xi] (pi, fi) ^
(m+u) I
\ ^m iXi;(j)M,fM) J
and consider the following identity for 2M x 2M determinant
(72)
^m V^m+M-3 ^m+M-1 Ym+M~2
^ Vm+M-l Ym+M-2 Ym-1
—(u-l)
Vm+M~2
—iv-l) —[v-l)
Vm+M-3 V'm+M-2
The identity is valid since this determinant can be transformed to
fli:^ ■■■ fm+M-3
—(I'-l) —(u)
fm+M-1 Vm+M-2
—{u~l)
H^m ^m+M-3 ^m+M-2
by adding and subtracting the columns and rows. Applying the Laplace expansion to the
l.h.s. of identity ( |7^ ) we obtain the bilinear relation for M x M determinants
'rm 'i^m+M-3 ^m+M-l V^m+M-2
(u-l) (^-l) _ _ (u-l) (u-l)
Ym~l '■fm 'Pm+M-3 'rm+M-2
rm Ym+M-3 V^m+Af-2 rm+M-2
(.-1) (,_1) _ _ (.-1) [u-l)
rm-1 rm rm+M-3 rm+M-l
^{u~l) . . .J-~^) ^(-1) ^(-1)
rm rm+M-3 rm+M-1 rm+M-2
(.-1) (,_1) _ _ (.-1) H
rm-1 rm rm+M-3 rm+M-2
0.
Using now the dispersion relations ( [70| ) we can identify
_(!/-!) — (j/-l) —{J^-I) — (I'-l)
V^m-1 V^m ■ ■ ■ '^m+M-3 Vm+M-2
—{v-1) — (I'-l) — ('^-l) — (i')
V^m ■ ■ 'Vm+M-3 fm+M~l fm+M-2
—(„-!) —(v-l) — (i^-1) —M
^m ' ' ' '^m+M-3 V^m+M-2 Vm+M-2
_(i/-l) — (jy-1) — (i^-1) —(1^-1)
V'm-l Vm ' ' ' '^m+M-3 Vm+M-l
—{u-l) —{v-1) —{u-l) —{u)
Vm-l Vm ' ' ' '^m+M-3 Vm+M-2
T,
{u-l)
m—1 ?
M
1 dr,
1 + K dxi
1
-H
VT
=r.
{^)
K
1 Srl'^-^)
'm— 1
\/l + K dxi
m—1
SO the bilinear identity for determinants becomes the bilinear identity for tau-functions
M
1 dr,
1 + K dxi
-H . J'^-'^)
T.
m—1
VT
=T,
M
K
Vl + fi: t^a;i
' m m—1 ^
20 S. PAKULIAK AND S. SERGEEV
which can be written in the form
^1 ['m-l " 'm J — K^ ^ i^) y-m ' m-1 ' m-1 'm
Last equahty coincide with (pSD at z/ = 0.
[1]
[2]
[3]
[4]
[5]
[6]
[7]
[8]
[9]
[10]
[11]
[12]
[13]
[14]
[15]
References
Pakuliak, S. and Sergeev, S. Relativistic Toda chain at root of unity. Preprint alin. SI/0101024
Pakuliak, S. and Sergeev, S. Relativistic Toda chain at root of unity II. Modified Q-operator. Preprint
ITEP-TH-39/01, lalin. SI/0107062
V. V. Bazhanov and N. Yu. Reshetikhin, "Chiral Potts model and discrete Sine-Gordon model at roots
of unity", unpublished (1995).
Bazhanov, V.V. and Stroganov, Yu.G. Chiral Potts model as a descendant of the six-vertex model, J.
Stat. Phys. (1990) 59, 799-817.
Faddeev, L.D. and Takhtajan, L.A. Hamiltonian methods in the theory of solitons, Berlin, Springer,
1987.
Newell, A.C. Solitons in mathematics and physics, Society for Industrial and Applied Mathematics,
1985.
Brusi, M. and Ragnisco, O. The periodic relativistic Toda lattice: direct and inverse problem. Inverse
Problem (1989) 5, 389-405; Lax representation and a complete integrability for the periodic relativistic
Toda lattice, Phys. Lett. (1989) A134, 365-370.
Kuznetsov, V.B. and Tsiganov, A.V. Separation of variables for the quantum relativistic Toda lattices.
Preprint hep-th/9402111
Suris, Yu.B. Discrete time generalized Toda lattices: complete integrability and relation with relativistic
Toda lattices, Phys. Lett. (1990) A145, 113-119.
Ruijsenaars, S. Relativistic Toda systems, Comm. Math. Phys. (1990) 133, 217-247.
Kharchev, S.; Mironov, A.; Zhedanov, A. Faces of Relativistic Toda Chain, Int. J. Mod. Phys. (1997)
A12, 2675-2724.
Ohta, Y., Kajiwara, K., Matsukidaira, J. and Satsuma, J. Casorati determinant solution for the rela-
tivistic Toda lattice equation, J. Math. Phys. (1993) 34, 5190-5204.
Jimbo, M. and Miwa, T. Monodromy preserving deformations of linear ordinary differential equations
with rational coefficients II, Physica D (1981) 2, 407; III, Physica D (1981) 4, 26.
Ueno, K. and Takasaki, T. Toda Lattice Hierarchy, Adv. Studies in Pure Math. (1984) 4, 1-95.
Date, D., Jimbo, M., Kashiwara, M. and Miwa, T. Transformation groups for soliton equations, in
Proceedings of RIMS sym. Nonlinear integrable systems - classical and quantum theory, eds., M. Jimbo
and T. Miwa, World Scientific, Singapore 1983, p. 39.
21
[16] Date, D., Jimbo, M. and Miwa, T. Solitons: differential equations, symmetries and infinite dimentional
algebras, Cambridge Tracts in Mathematics, 135, University Press, Cambridge, 2000.
[17] Hirota, R. and Satsuma, J. Theor. Phys. Suppl. (1976) 59, 64.
BoGOLiuBOv Laboratory of Theoretical Physics, Joint Institute for Nuclear Research,
DuBNA 141980, Russia and Max-Planck Institut fijr Mathematik, Vivatsgasse 7, 53111 Bonn,
Germany
E-mail address: pakuliak@tlisuiil.jinr.ru
BoGOLiuBOv Laboratory of Theoretical Physics, Joint Institute for Nuclear Research,
DuBNA 141980, Russia
E-mail address: sergeev@thsunl.jinr.ru